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Abstract 

Lepton masses and mixing angles via localization of 5D fields in the bulk are revisited in the context of 
Randall-Sundrum models. The Higgs is assumed to be localized on the IR brane. Three cases for neutrino 
masses are considered: (a) The higher dimensional LH.LH operator (b) Dirac masses (c) Type I see-saw 
with bulk Majorana mass terms. Neutrino masses and mixing as well as charged lepton masses are fit in 
the first two cases using y 2 minimisation for the bulk mass parameters, while varying the 0(1) Yukawa 
couplings between 0.1 and 4. Lepton flavour violation is studied for all the three cases. It is shown that 
large negative bulk mass parameters are required for the right handed fields to fit the data in the LH LH 
case. This case is characterized by a very large Kaluza-Klein (KK) spectrum and relatively weak flavour 
violating constraints at leading order. The zero modes for the charged singlets are composite in this case 
and their corresponding effective 4-D Yukawa couplings to the KK modes could be large. For the Dirac 
case, good fits can be obtained for the bulk mass parameters, Cj, lying between and 1. However, most of 
the 'best fit regions' are ruled out from flavour violating constraints. In the bulk Majorana terms case, we 
have solved the profile equations numerically. We give example points for inverted hierarchy and normal 
hierarchy of neutrino masses. Lepton flavor violating rates are large for these points. We then discuss 
various minimal flavor violation (MFV) schemes for Dirac and bulk Majorana cases. In the Dirac case with 
MFV hypothesis, it is possible to simultaneously fit leptonic masses and mixing angles and alleviate lepton 
flavor violating constraints for Kaluza-Klein modes with masses of around 3 TeV. Similar examples are also 
provided in the Majorana case. 
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I. INTRODUCTION 



One of the most interesting solutions of the hierarchy problem is the Randall-Sundrum model 
[T] which proposes a warped extra space dimension compactified on an S\jZ% orbifold. Two branes 
representing the UV and the IR scales are located at the two end points of the orbifold. In the 
simplest models, the Standard Model matter and gauge fields are localized on the IR brane along 
with the Higgs field. Massive Planck scale modes are exponentially suppressed at the IR brane, due 
to the warped bulk geometry, caused by the presence of a large negative cosmological constant 1 . 
Variations of this set up have been considered in several different contexts 2 . 

For example, introducing gauge fields in the bulk facilitates unification of couplings |4J. But 
this leads to large corrections to the electroweak precision observables and places a lower bound 
on the mass of the lightest gauge Kaluza-Klein (KK) mode to be around 25 TeV. This is because 
the coupling of brane localized fermions to the gauge KK states is enhanced by a factor ~ 8.5 
compared to the SM coupling [5HZ] ■ A similar study in terms of oblique parameters was reported in 
[8, 9j. Boundary kinetic terms for the gauge fields can lower the bound [101111], but this might spoil 
the unification. Alternatively, allowing the fermions to propagate in the bulk eases the constraint 
of 25 TeV on the lightest KK mode, to about 10 TeV [12]. Having a bulk Higgs further eases 
the bound [Jj. On the other hand, scenarios with extended particle content and a bulk custodial 
symmetry with a brane localized Higgs boson were found to lower the bounds on the KK gauge 
boson mass to ~ 3 TeV [13]. In p3] the authors explored a mixed scenario where part of fermions, 
the third generation quarks are localized on the IR brane. It was shown that such a scenario would 
soften the corrections to the p parameter. Finally modifying the RS metric near the IR boundary 
can also help in reduction of the strong electroweak precision constraints [151 ES] • 

Allowing fermions to propagate in the bulk has interesting implications for flavor physics. The 
bulk profiles of the fermion fields are determined by their bulk masses in a manner similar to 
Arkani-hamed and Schmaltz mechanism in ADD models [T7]. In the RS model, however, the 
warped geometry facilitates the so-called 'automatic' localization of fermions [6j. The profiles are 
also no longer gaussian, but are exponentially suppressed. It has been proposed that RS could 
be a theory of flavour, where the fermion mass hierarchy can be explained in terms of a few 

1 The RS metric is given by 

ds 2 — e~ 2,y( - y - ^rj^dx^ 1 dx v — dy 2 , 

where cr(y) — k\y\. For recent reviews on RS models, please see [2j. 

2 The phenomenology of RS models has been extensively studied. A recent review on collider phenomenology 
concentrating on LHC can be found in [3J. 
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0(1) parameters. This is analogous to the popular Froggatt-Nielsen (FN) models [181 EH] m f° ur 
dimensions. While in the FN model, it is the gauge and the heavy fermion sector which determine 
the hierarchies in the Yukawa couplings, in the RS case, it is the geometry of the bulk. The role 
of the FN charges can be played by the five dimensional Dirac masses for the bulk fermions. The 
expectation is that by taking 0(1) bulk mass parameters as well as Yukawa couplings, one would 
be able to explain the large hierarchies in the quark and leptonic mass spectrum. While this is true 
in general for quarks and charged lepton masses, as we will see subsequently, in case of neutrino 
masses, the situation is a bit more involved. 

Flavour violation in the hadronic sector has been explored by various authors |2UH22| . a recent 
comprehensive analysis can be found in [23(. 124] . In the present work, we are interested in studying 
neutrino masses and mixing angles within the RS context. One method of generating neutrino 
masses in the RS model would be to allow only the right handed neutrino to propagate in the 
bulk, while the SM particles are confined to the IR brane. This leads to a higher dimensional 
seesaw mechanism [25]. However, unlike the case of ADD models, here only the lightest KK 
modes participate in the seesaw mechanism. Furthermore, lepton flavour violating decay rates are 
extremely large in this case pushing the lightest KK mode to be heavier than m^k > 25 TeV [26] . 
Neutrino mass models have also been explored in the alternative scheme where all the fermionic 
fields are allowed to propagate in the bulk. In the present work, we will concentrate on this set up 
and study the neutrino mass phenomenology and lepton flavor violation [25 - 33] . We have assumed 
Higgs to be localized on the IR brane. Fermion mass fits in scenarios with Higgs also propagating 
in the bulk have been considered in |21| IM] . 

In this RS set-up (fermions in the bulk, Higgs localized on IR brane ) neutrino mass models 
can be divided broadly into Dirac mass models or Majorana mass models. In the case of Majorana 
fermions, the number of possibilities is more than one. In the present work we discuss three cases 
in detail (a) The higher dimensional LH LH operator (b) the Dirac neutrino case and finally (c) 
Majorana neutrinos with bulk seesaw terms. In these models, typically two sets of parameters 
determine the charged lepton masses and neutrino masses and mixing angles. These are the afore 
mentioned set of bulk Dirac masses for the fermions and then the 0(1) parameters containing the 
Yukawa couplings. In each of these cases, we have numerically minimized a x 2 function containing 
the model parameters and the leptonic masses and mixing data, to determine the 'best fit' regions 
of the parameter space. The Yukawa couplings are varied from 0.1 to 4 whereas the ranges for the 
bulk parameters are judiciously chosen to be as wide as possible. 

We found that in the (a) higher dimensional LHLH operator case, the bulk mass parameters 
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of the charged singlets are required to be negative and extremely large. This gets reflected into 
an extremely hierarchal Kaluza-Klein mass spectrum of the first KK states of the SM fermions. 
In fact, the best fit regions are those with Standard Model charged singlets being completely 
composite 3 . On the other hand, if one considers Dirac neutrinos, it is quite possible to fit the data 
naturally with the bulk Dirac masses within reasonable ranges without any large hierarchies. Both 
hierarchal and inverse hierarchal neutrino mass schemes can be fit in this case though it is much 
more difficult to find regions which satisfy inverse hierarchal neutrino mass relations compared to 
normal hierarchy. The bulk equations of motion in the presence of a Majorana mass term are 
coupled and more complicated than the Dirac or LHLH case. We have solved them numerically 
and given example points where data can be fit easily either the inverted or the normal hierarchy 
scheme. We have not conducted an extended numerical scan of the parameter space for the bulk 
Majorana case. 

Fitting neutrino masses in any of the above models in RS set up potentially leads to large 
lepton flavor violation. A detailed analysis was presented in [30], where the authors discussed 
the implications of flavor physics in the lepton sector with both the brane localized and the bulk 
Higgs. Neutrinos were assumed to be of Dirac nature. They observed that with a bulk Higgs, the 
branching fraction for the process fj, — > &y requires a KK mass scale of around ~ 20 TeV to keep 
it below the present experimental limits. Similar comments were made in |33j on how the higher 
dimensional operator case is not conductive for suppressing process like [i —> eee, especially when 
the KK mass is low. Higgs was allowed to propagate in the bulk in this work. In the present work, 
we revisited the flavor constraints for all the three cases, concentrating on the best fit regions in 
the LHLH and the Dirac case. For the LHLH case, the couplings of SM fermions to KK gauge 
bosons are universal in the best fit region, leading to no apparent constraint, at least at the leading 
order from the tree level flavor violating decays. However, there are large Yukawa couplings in this 
model which make it unattractive from perturbation theory point of view. The best fit region of 
the Dirac case is strongly constrained from tree level decays as well as loop induced decays like 
H — >• e + 7. In the brane localized Higgs scenario we are considering here, the limits from dipole 
processes are cut-off dependent. But, for cut-off values close to the first KK mass scale, the limits 
are comparably much stronger. For the bulk Majorana case too, the points we have considered 
display strong constraints from leptonic flavor violation and are ruled out. One would thus need 
ways to circumvent these strong limits from lepton flavor violation. 

3 This interpretation is based on the AdS/CFT correspondence. 
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We explored Minimal Flavour Violation (MFV) ansatz implemented in the RS scenario to evade 
the flavour constraints in the Dirac and Majorana cases |31| 135] . We provide example symmetry 
groups where the flavor violating constraints can be removed for both the Dirac and the Majorana 
cases. 

The paper is organized as follows. In section ([II]), we discuss lepton mass fits in three models 
of neutrino mass generation, the higher dimensional LHLH operator, the Dirac case and the bulk 



Majorana mass terms case spread over three subsections. In section (III), we discuss the lepton 



flavor violating constraints for the three cases of neutrino masses. In section (IV) we discuss the 
minimal flavor violating schemes for the Dirac and Majorana cases and show example points where 
flavor violating constraints are alleviated. We close with a summary and outlook in the final section 

m 



II. LEPTON MASS FITS 



The observed neutrino and charged lepton data is fit to the set of theory parameters which 
determine the charged lepton and neutrino mass matrices through a x 2 minimization. Thus the 
observables correspond to three charged lepton masses, three mixing angles and two (neutrino) 
mass squared differences, while, the bulk mass parameters and Yukawa couplings form the set of 
theory parameters. The number of theory parameters varies from model to model, as discussed in 
the following sub-sections. We have chosen the following central values for the observables [36, 37 : 



TABLE I: Experimental Data 



masses 
(MeV ) 


mass-squared 
(eV 2 ) 


mixing angles 


m - fl M +0.0000007 
'">e — u - dl -0. 0000007 

777 - lf)^ fi + °- 000003 
M — ruo.u_ 000003 

777 - 1 776+0-00016 
"'t — x ' ' °-0. 00016 


Am 2 12 = 7.59±g;|S x 1(T 5 
Am§ 3 = 2A3t° il x 10 ~ 3 


012 = 0.59i°;°i5 

023 = o.79±g:i| 

13 = 0.1541°;°^ 



We use the standard x 2 definition for N observables given by 

N / exp theory \ ^ 
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where, y i eory is the value of the i observable predicted by the model and y\ xv is its corresponding 
experimental number measured with a uncertainty of <7j. Since, the values of the charged lepton 
are measured to a very high accuracy, it is difficult to fit masses to such high accuracy. Thus, we 
incorporate up to ~ 1.5% errors in the masses of charged leptons 4 . The x 2 relevant to our study is 



2= (^-0.59) 2 (e atm -0.79) 2 (fl 13 -0-154) 2 (Am 2 sol - 7.59 x 10~ 23 ) 2 
X (0.02) 2 (0.12) 2 (0.02) 2 (0.2 x lO" 23 ) 2 

(Am 2 atm - 2.43 x 1(T 21 ) 2 (m e - 0.00051) 2 (m M - 0.1056) 2 (m T - 1.77) 2 
+ (0.2 x 10- 21 ) 2 + (0.00001) 2 + (0.0001) 2 + (O02) 2 () 

As mentioned above, the fermion masses (and mass squared differences) and mixing angles 
appearing in Eq.Q are functions of bulk parameters. The minimization was performed using 
MINUIT |40j . For a given scan, MINUIT looks for a local minima for the x 2 around a certain 
input guess value of the bulk masses and Yukawa parameters. This scan is repeated by randomly 
varying the guess values and in the process of looking for a global minima. 



A. The LHLH operator 

In the absence of detailed specification of the mechanism which generates neutrino masses, one 
can always write an effective higher dimensional operator at the weak scale to account for non-zero 
neutrino masses. In the Standard Model, this operator is simply the (LHLH)/ A operator , where 
A is the high scale at which neutrino masses are generated. In the Randall Sundrum model one 
can write a similar operator for non-zero neutrino masses. The model has been earlier studied in 
|21| [28] . The 5D action for the RS model with the Higgs localized on the IR brane is given by 

g ( 1(0 - m L )L + E(0 - m E )E ) 

g {^LHLH + Y E LEH) 5(y - ttR) (3) 

where ~ 2.2 x 10 18 GeV is the fundamental five dimensional reduced Planck scale and 

D M = d M + n M + l -fr a W^(x, y) + |-Qy5 M (x, y) (4) 

with Ojvf = (—k/2e~ ky j^j 5 ,0) being the spin connection and Qy is the hypercharge. M is the 
five dimensional Lorentz index. R is the compactification radius and k and Ye are the coupling of 

4 This approach is very similar to fermion mass fitting in GUT theories. See for example, |38l I39| . 



S = Skin + 5Vuk 

S kin = f*xfdy 
SVuk = d 4 x dy 
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the neutrino mass operator and the Yukawa coupling for the charged leptons respectively. They 
are three dimensional matrices in flavour space and we have suppressed the generation indices in 
writing the above equation. L and E are the 5D fermionic fields which transform as doublets and 
singlets respectively under the Standard Model SU(2) W gauge group with the covariant derivative 
given by Eq.Q acting accordingly, mi, and the are five dimensional Dirac masses of the L and 
E fields. As we will see below, after Kaluza-Klein decomposition, these masses determine the 
profiles of the zero and higher KK modes in the extra dimension. Since the effective operator is 
suppressed by the 5D Planck mass, one can imagine that the neutrino masses are as a result of 
some fundamental lepton number violation beyond the 5D Planck scale. 




FIG. 1: Regions in Cj for the LHLH case which give best fit to lepton masses and mixing. The graphs in the 
upper row shows the region of parameter space for the bulk masses for doublets which fits small neutrino 
masses. Neutrino masses are assumed to have normal hierarchy in this analysis. The graphs in the lower 
row shows the region for the bulk masses for the charged singlets c^. We have used log scale for cg^ 



The left and right components of the L and E fields have different Z2 properties. These are 
chosen such that the Z2 even zero modes correspond to the SM fields. We assign the following Z2 
parity for the L^ r and the E^ r fields, where the subscript (7,r) correspond to the left and right 
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handed components of L and E 5 . 

Z 2 (y)Li(x,y) ->• Li(x,y) , Z 2 (y)L r (x,y) ->■ -L r (x,y) 
Z 2 (y)E r (x,y) -)• E r (x,y) , Z 2 {y)E l (x,y) -> -E t (x,y), 

where Z 2 (y) ■ y —> —y. The 5D fields can be expanded in terms of the KK modes, with the 
expansion given by (20J [25] 

^(x,y) = X:^=e 2 ^)^(x)/f ) ( y ) ; L r (x,y) = X:^ D e 2 ^)4n)( x ) x (™)( 2/ ) 

n=0 n=0 
oo 1 oo 1 

= J3 ~i=e^E^{x)ft\y) ; = E -i=e^E { r\x) X ^\y) (5) 

n=0 n=0 

where the exponential factor is chosen such that the fields are canonically normalized. The profiles 
fL,E and Xl,e are determined by : 

(dy + c L a')f^ E (y) = m^e^ X ^ E (y) 
{-d y + c L a') X ^ E (y) = m^e^fi%(y) (6) 

where the 5D masses m,L,E are written in terms of the fundamental scale as mjj^; = cl^e^' and 
a' = dy<j = k. The following orthonormality conditions are used for the profiles jh,E and xl,e to 
arrive at Eq.Q 

1 r * d y e^ E {y)xtky) = -7^ r d v e°&(y)ftMy) = * nm (?) 



wR 



The above equations decouple for the zero mode solutions where = 0. The solution for the 
Z 2 even part, /z,(y) is given as 



f?\y) = N (c L )e-^'y ; N (c L ) = j! 2 J^t_ 1 (8) 

iVo being the normalization constant. The solution is the same for profile of E, with 
cl replaced by eg. The bulk wave functions are exponentials which peak towards the UV (IR) 
for c > 1/2 (c < 1/2) as can be seen from Eq.Q. Typically, particles lighter in mass like leptons 
require c > 1/2 whereas heavier particles like top quark is localized much closer to the IR brane with 
c < 1/2. For the charged leptons and the neutrino masses one would expect all the corresponding 
Cj to be > 1/2. The KK expansions ^ are put into the Yukawa part of the action Eq.Q leading 

5 The 75 required to define the left and right components remains the same as the four dimensional case. 
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to 



S- 



Yuk 



d A X 



( Y E L^(x)fl n \y)E^{x)fP{y)e kR -H 



+ 



(9) 



where we have used — > e H to canonically normalize the Higgs field and suppressed the 
subscripts (I, r) for the Z2 even fields. The odd fields are neglected as they are removed from the 
boundary as a consequence of the Z2 symmetry. The charged lepton mass matrix and the neutrino 
mass matrix are determined when the zero modes of the fields are taken. The charged lepton mass 
matrix, corresponding to the L^'E^H operator in the action is given by 

„3 



Y E 



Ye 
Rtt 



No(c L )N (c E ) e 



(1— Cf,— CE)kRw 



(10) 




FIG. 2: The distribution of electron Yukawa couplings (Ye) w hich give a 'good fit' to the charged fermion 
mass data in the LH LH operator case. Neutrinos are assumed to follow normal hierarchy in this analysis. 
The binning is done with an interval of 0.2 



where the matrix Ye can be considered equivalent to the 4D dimensionless Yukawa couplings. 
The neutrino mass matrix defined as the co-efficient of the L^L^HH operator in the action, is 
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given as 

M »v ~ 3 2A( 5 ) 1 m kk y A( 5 ) 

= «y e 2kR *f Li (nR)f L] (7rR) = g ^(czjA^.) e (^-^)«* (n) 

where i, j are generation indices and Mkk is the typical mass of higher KK fermions. The correc- 
tions are from higher order KK modes and can be neglected. Before fitting the mass matrices, we 
introduce new 0(1) Yukawa parameters entering the mass matrices, which are defined as 

Y E = 2kY E ; k =2kn (12) 

In terms of these new Yukawa parameters, the mass matrices are explicitly given as 



(M(°>°))y = —(Yi;)^ 1 -™^ 



(0.5 -c Li ) / (0.5 -c E . 



V2 J V e^~ 2cL > kR -l\ e ^-^E j )M_ 1 



f y(o,o)v. _ v2 , j, {2-c L -c L .)kR« (0-5 -c Li ) / (0-5 -c Lj .) 

The matrices are diagonalised as J7| L 7Wi°' ^C/ e i{ = Diag[{m e , m^, m T }] and [/^.M^ ° £/J = 
Diagfjm^, m U2 ,m Ki }] and Upmns = utu e L- The eigenvalues of the charged lepton mass ma- 
trix and the mass squared differences of the neutrino mass matrix and the Upmns mixing angles 
are fit to the data as per Table [TJ In this case, there are three cl and three c E - and fifteen Yukawa 
parameters fitting three charged lepton masses, three angles and two mass squared differences. 
Given the dependence of the leptonic mass matrices on the Yukawa parameters, we have chosen 
them strictly to be of 0(1) nature. By this we mean, they are varied roughly between -4 and 4. 
Furthermore, in order to avoid regions where the Yukawa parameters are unnaturally close to zero, 
we put a lower bound on the Yukawas such that |Y| lies between ~ 0.08 and 4. 

Since the charged leptons and neutrinos have relatively light mass spectrum compared to heavy 
quarks, one would have expected that varying cl and c E between 1/2 and 1 would be sufficient 
to fit the data. However, in the present context such values for c E will not satisfy the data. This 
is because the neutrino mass matrix depends only on c Ei and requiring the neutrino masses to be 
of the 0(lO _1 )eV automatically sets CL i to be around 0.9, close to the UV brane. The charged 
lepton mass matrix, which in turn is determined by both c^. and c Ei should off-set the effect of c Ei 
and increase the effective 4D Yukawa coupling by pushing it towards the IR brane. This can only 
be achieved by taking large and negative values 6 of the c Ei . The range for the scan of the cl : e 



6 One way to avoid large negative c parameters would be to consider very large 0(1) Yukawa parameters. The 
required Yukawa couplings are in the range ~ O(10 3 — 10 4 ) to make any connection with data. 
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-n-r-^i-r-r-mTTl 



-htttttt-tttt-H 



0^ 



rrnrrrrr^TT-rr^rrr 




FIG. 3: The distribution of neutrino Yukawa couplings («/) which give a 'good fit' to the fermion mass data 
in the LH LH operator case. Neutrinos are assumed to follow normal hierarchy in this analysis. The binning 
is done with an interval of 0.2. 



has been judiciously chosen between 0.82 and 1.0 for bulk doublets and —5 x 10 7 < ce 1 < —0.2, 
— 10 8 < ce 2 < —8000 and — 10 9 < ce 3 < —9000 for first, second and third generation charged 
singlets respectively. A larger democratic range does not change the results significantly. 

All the parameters, the fifteen Yukawa couplings and the six cl,e parameters are varied so as 
to minimize the function in Eq.Q. The points which give a x 2 between 1 and 8 are considered 
to give a 'good fit' to the data. In Fig.jlj we present the regions in cl 1 2 3 and ce 1 2 3 which have 
minimum \ 2 assuming normal hierarchy for neutrino masses. It is important to remember that 
Yukawa couplings are also varied in obtaining this range in the cl e parameter space. From the 
figures we see that the strong constraint of neutrino masses limits the c/^ to be within a limited 
range. On the other hand, ce seem to have much larger ranges spanning orders of magnitudes. 
In particular, ce 1 is virtually unconstrained from 0{— 1) to 0{ — 10 6 ). This is an artifact of the 
unconstrained lightest neutrino mass, m ui . ce 2 and ce 3 have lesser freedom as they are constrained 
by the mass squared differences. The allowed ranges in the cl,e which satisfy the minimum \ 2 
requirement are summarized in Table [11} 

It would be interesting to see distribution of the Yukawa couplings Y' E and k! for the 'best fit' 
regions of the parameter space. The distributions are presented in Figs. [2] and [3j. For most of 
the YL parameters, there is peaking at the two ends of the range chosen, around 0.2 and 3.8. The 
exception is the lower 2x2 block of the Yukawa matrix, for which there seems to be a flatter profile 
for the upper row parameters {Y'e)^ ana ^ C^e)23 and a progressively increasing distribution for the 
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TABLE II: Allowed range for the bulk parameters with minimum \ 2 ■ Neutrino masses have normal hierarchy. 
Range of first KK scale of the doublet (singlet) M^\m^) corresponding to the bulk mass parameter is also 
give. 



parameter 


range 


ranj 


;e of M ( L 1] (TeV) 


parameter 


range 


range of M^ 1} (TeV) 


CLi 


0.87-0.995 




1.49-1.59 


CEi 


-10.0 to -5.0 x 10 6 


7.9-3.9 x 10 6 


CL 2 


0.86-0.98 




1.48-1.58 


ce 2 


-1.0 x 10 4 to -1.2 x 10 8 


7.9 x 10 3 -9.5 x 10 7 


CL 3 


0.84-0.92 




1.47-1.53 


ce 3 


-7.0 x 10 5 to -1 x 10 9 


5.5 x 10 5 7.9 x 10 8 



second row parameters. 

For almost all of the Y' E parameters, peaking seems to be happening at high values ~ 3.8, 
except for {Y' E )22- There are also second peaks at very low values ~ 0.2 for some of the parameters. 
Distributions in k' on other hand, show peak at very large value ~ 3.8 for the first two generation 
couplings and very low values ~ 0.4 for k' 33 and k' 23 . With the exception of peaks, there is an 
underlying though highly subdued, 'anarchical' nature in the distribution of Y' E Yukawa couplings 7 . 
Thus, for a given choice of 0(1) Yukawa couplings within our chosen range (-4 to 4), it seems to 
be possible to find c values which can fit the data well 8 . From the allowed parameter space, we 
have randomly chosen two sample points, which we call Point A and Point B, and we provided 



the corresponding observables in Table III The corresponding Yukawa couplings are given in Eqs. 



(14) and (15) 



Yukawa coupling matrices for Point A: 



0.5023 1.9546 3.9730 
3.2482 2.9629 2.7742 
2.6865 2.0383 1.2369 



Yukawa coupling matrices for Point B 



Y'e 



3.0571 0.6316 0.8978 
1.4085 0.9952 3.5597 
0.7971 0.9579 0.5539 



3.8933 3.9717 3.9818 
3.9717 -2.6660 -1.1409 
3.9818 -1.1409 1.4597 



0.2315 -3.8320 0.3490 
-3.8320 -0.6632 -1.1287 
0.3490 -1.1287 0.0802 



(14) 



(15) 



7 Anarchy in the Yukawa distributions does not necessarily mean anarchical structure in the mass matrix. 

8 Increasing the scan range for the 0(1) Yukawa couplings from -10 to 10 does not change the gross features of the 
distributions much. For example, Y' E are peaked near the end points, showing that the lepton masses in this case 
prefer large or small Yukawa couplings. The re' distribution has the same features scaled now to to to 10 from 
to 4. The ranges of the cl,e do not change significantly. 
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TABLE III: Sample points with corresponding fits of observables for Normal Hierarchy in LHLH case with 
0(1) Yukawas. The masses are in GeV 



Point 


A 


B 


x 2 


2.07 


5.5 




0.9755 


0.903 


CL 2 


0.9162 


0.93 


cl 3 


0.87 


0.8443 


Ce 1 


-692416.99 


-17.35 


ce 2 


-2647794.18 


-946125.13 


Ce 3 


-80717122.21 


-47941542.53 


m e 


5.07 x I0~ 4 


5.08 x 10- 4 


"V 


0.1056 


0.1056 


m T 


1.767 


1.771 


#12 


0.58 


0.589 


#23 


0.68 


0.743 


#13 


.168 


0.163 




7.49 x 10~ 23 


7.48 x 10~ 23 




2.47 x 10~ 21 


1.99 x 10~ 21 



In Appendix [A] we have presented our results assuming neutrinos have an inverse hierarchical mass 
ordering. We find very few points which satisfy the data in this case. This is because inverted 
hierarchical spectrum requires two masses at the atmospheric neutrino scale with their mass dif- 
ference satisfying Am 2 ol . Thus the results are very sensitive to the 0(1) Yukawa parameters. For 
a fixed Yukawa, however it is easy to find points. More discussion is present in Appendix [Aj 

The analysis presented so far has been purely phenomenological. Let us digress from the fermion 
fits for a moment to discuss about the large negative c parameters. Such large negative values for 
the bulk mass parameters are in conflict with the 5D cutoff scale k. We have neglected this conflict 
in fitting the data where we have considered them to be purely phenomenological parameters which 
can take any value 9 . In terms of the bulk wave- functions the large negative c values would mean 
that the zero mode wave- function 3> 1, which is not the case when we choose the c parameters 

9 We prefer to keep the k (and also the radius R) value fixed by noting that only the charged singlets required large 
negative c values. In case we shift the 5D cut-off scale to \c\k keeping k fixed, the corresponding IR would shift to 
cAir, thus spoiling the solution to the hierarchy problem in this scenario. 
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between and 1. 

It is preferable to understand the large negative c values in terms of localization on the IR 
brane. The limit c — > — oo corresponds to the case where the fermions are completely localized 
on the IR singular point [41J. In the limit c — )• — oo, /(°) — > oo indicating full overlap of the bulk 
wave- function with the brane. The value of the c parameters also affects the masses of the KK 
modes. These masses are determined from Eqs.Q by considering m n ^ and choosing appropriate 
boundary conditions for the 5D fields. The resultant differential equation has solution in terms 
of Bessel's function which describe bulk wave-functions of the KK modes whereas the masses are 
given in terms of the zeros of the Bessel function [42] . The order of the Bessel function is roughly 
given by |c| for large values of c. In the asymptotic limit the first KK mode has mass ~ \c\ke~ kRn . 
Thus we see that the phenomenologically relevant first KK mode mass also grows as ~ cAjr, where 
Ajr ~ TeV, the IR cutoff. The masses of the first KK modes are presented in the Table [il]. The 
bulk wave- function of the KK mode tend to zero as |c| — > oo. 

One might wonder if such large negative values of the parameters would have some im- 
plications in terms of the AdS/CFT correspondence [2} I43j. The CFT interpretation for the bulk 
scalars has been studied in [21 [44] and for bulk fermions in |45j . The best fit cl,e parameters 
of LLHH case given in Table [II] leads to an unusual situation where the left handed leptons are 
almost completely elementary while the right handed singlets are completely composite. This can 
be easily verified using the 'holographic basis' of [46]. The composite component of the cl is pro- 
portional to e -( c L-o.5)kRTT ^ w j 1 j c ] 1 g oes to when cl —> 0.99. Thus, the zero modes for the doublets 
are elementary. For the ce fields however, the elementary component for the zero mode is given 
as yj (ce — 0.5) (ce + i.5)e~l 1,5+C£; l fc ' R7r . Thus we see that the zero mode for the charged singlets 
have a vanishing elementary component and are completely composite fields. The effective 4-D 
Yukawa coupling of the zero mode to the KK modes, is given as Y' E \J (0.5 — ce). A problematic 
feature of these models is that this coupling enters the non-perturbative regime for ce large and 
negative. This non perturbative coupling appears for all including the first KK mode, which is 
phenomenologically relevant. This, non-perturbative feature is restricted to the Yukawa coupling. 
The gauge coupling on the other hand do not face this problem. In fact as we shall see later (Sec- 
tion V, Figure [l3]), the coupling strength of the zero mode fermions to gauge KK modes quickly 
approaches the coupling of the brane localized fermions to gauge bosons for relatively moderate 
values of |c| parameter. 
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FIG. 4: The figures above correspond to the case in which neutrinos are of Dirac type. The points in the 
above figures correspond to a x 2 between 1 and 8. The plot represents the parameter space for the bulk 
masses of the doublets. This case corresponds to the normal hierarchial case. 



B. Dirac Neutrinos 

Dirac neutrino mass models in the RS setting have been extensively studied in the literature 
|30| . In [33], the authors talked about the difficulty of fitting neutrino masses and mixing angles 
in the same scenario as quarks. Their argument drew inspiration from the fact that neutrino 
mixing angles are anarchic in nature. To address this issue they had a bulk Higgs, with the 
profile 'sufficiently peaked' near the IR brane and introduced a 'switching behaviour' to fit the 
both charged fermion and the neutrino masses and mixing angles. We, on other hand, approach 
this problem in the same way as we have done in the LHLH case of the previous section. We 
look for regions in the parameter space of the bulk masses which give 'good' fits for a reasonable 
choice of 0(1) Yukawa couplings. The particle spectrum of the Standard Model is extended by 
adding singlet right handed neutrino. Global lepton number is assumed to be conserved. It can be 
violated by quantum gravity effects which manifest at the 5D Planck scale. However, for most of 
the present analysis, we require lepton number violation present to be highly suppressed. 

The bulk and Yukawa actions in Eq.([3j) now take the form: 

Skin = J d 4 xj dy^j ( L{0 - m L )L + E{0 - m E )E + N{0 - m N )N ) 

S yuk = Jd^xjdy^i Y N LNH + Y E LEH ) 5(y - nR), (16) 
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FIG. 5: The plot represents the parameter space for the bulk masses of charged singlets. 



where N stands for the 5D right handed neutrino fields. The rest of the parameters carry the 
same meaning as in the previous section. The components of the N field are assigned the same Z<i 
properties as the E field. We expand the N fields as 

oo 1 oo 1 

N r (x,y) = V -=e 2 ^N^(x)fft\y) ; Nfay) = V N[ n \x)xP (y) (17) 

" ''Kit n V Kit 



n=0 v n=0 

Using Eq.(17) and Eq.([5]) one can derive the equations of motion and solutions similar to Eq.Q 
for the profiles of N fields. Substituting them, the zero mode mass matrices for the charged lepton 
and neutrinos take the form: 



M [0fi) = Vy E ;Y E = ^ N (c L )N (c E ) e (^-c E )km 



^(0,0) 



(18) 



where we have neglected corrections from higher KK modes. As before, we perform a scan over 
the parameter space of the bulk fermion masses and order one Yukawa parameters to minimize the 
% 2 in Eq.|2| for the masses and mixing angles. To specify the parameters which are scanned, it is 



useful to look at the explicit form of the mass matrices equivalent to those of Eq.(13): 



(M°' 0) k 



(M°' 0) ) 



V 

^2 



(Y E )ije 



(l-c Li -c E )kRn / (0-5 ~ C Li 



(0.5 -c Ej ) 



e (l-2c Li )irkR _ 1 V e (l-2c E .)KkR 



V (V' \ ^-CLi-CNAkRir 



(0.5 - c u 



(0.5 -c N .) 



e (l-2c Lt )7vkR _ 1 V (l~2c Nj )n 



kR 



(19) 



where YL N = 2kY E ,N- Each of the Cj parameters (i = {L, N, E}) which are three in number are 
varied along with eighteen 0(1) Yukawa parameters, i.e, a total of 27 parameters are varied to fit 
the data and minimize the x 2 - The c parameters are varied as follows: The doublets (c^J and the 
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the charged singlets are varied between 0.02 and 1, while the neutral singlets are varied between 
between 1 and 1.9. The order one Yukawa couplings, Y EN , are varied randomly between -4 and 
4 with a lower bound \Y\ > 0.08. We consider all the regions of the q parameter space where the 
X 2 is between 1 and 8 as a 'good' fit region. In Figs. 4|5|6 we present regions in the c, parameter 
space which give 'good' fit to the leptonic mass and mixing angles. A summary of these regions is 



presented in Table (IV) 



TABLE IV: Allowed ranges of bulk parameters with normal hierarchy of neutrino masses. The range of first 
KK scale corresponding to the range of c values is also given. 



parameter 


range 


m£ ] TcV 


parameter 


range 


il4 1} TeV 


parameter 


range 


Ml 1] TeV 


CLi 


0.05-0.76 


0.839-1.4 


ce ± 


0.2-0.88 


0.959-1.5 


CjVi 


1.1-1.9 


1.67-2.31 


CL 2 


0.05-0.72 


0.839-1.37 


ce 2 


0.05-0.73 


0.839-1.38 


Cn 2 


1.1-1.9 


1.67-2.31 


°L 3 


0.05-0.64 


0.839-1.31 


ce 3 


0.05-0.64 


0.839-1.31 


cn 3 


1.1-1.9 


1.67-2.31 
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FIG. 6: The plot represents the parameter space for the bulk masses of the neutrino singlets. 



The Dirac neutrino mass matrix in the RS model seems to fit the data more naturally compared 
to the LHLH discussed in the previous subsection. A large section of the points fall in the regime 
Cj > 1/2 indicating that they are localized closer to the UV brane. The distributions of the Yukawa 
couplings in the 'good fit' region, presented in Figs.(7j8) show that most of them peak in the last 
bins for all the Yukawas at (3.8 — 4.0). A secondary peak can also been seen at (0.2 — 0.4) bin 
for some of the Y^ parameters. Electron Yukawa couplings on the other hand do not seem to 
show any such secondary peak. In this case too the distribution of the 0(1) Yukawa couplings 
displays an underlying anarchic nature especially for the Y' E . This will prove useful in our analysis 
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of Minimal Flavour violation where the 0(1) Yukawa couplings and the bulk mass matrices need 
to be simultaneously diagonalizable. In Table (|V|), we presented two sample points. Point A has 
all the Cj > 1/2 where as Point B has ce 2 ,ce 3 < 1/2. The corresponding Yukawa couplings are 
given m Eqs. pOpIj ). 

As before we use the holographic basis to comment on the partial compositeness of the bulk 
fermions. The zero modes of singlet right handed neutrinos are dominantly elementary, with almost 
zero component of compositeness. The composite component for the zero modes of the doublets 
and the charged singlets becomes smaller as the corresponding c values becomes greater than 0.5. 
Essentially they have partially composite nature. 



TABLE V: Sample points with corresponding fits of observables for Normal Hierarchy in Dirac case with 
0(1) Yukawas. The masses are in GeV 



Parameter 


Point A 


Point B 


x 2 


0.28 


0.39 


CLi 


0.6263 


0.7166 


C L 2 


0.5932 


0.6382 


CL 3 


0.5293 


0.6126 


CE X 


0.6704 


0.5911 


CE 2 


0.5541 


0.1939 


CE 3 


0.5131 


0.2647 


CNi 


1.2233 


1.2791 


Cn 2 


1.2692 


1.1215 


c Nl 


1.2948 


1.2343 


m e 


5.09 x 10~ 4 


5.09 x 10~ 4 


«v 


0.1055 


0.1055 


m T 


1.77 


1.77 


#12 


0.59 


0.589 


#23 


0.80 


0.792 


#13 


0.153 


0.153 


5m tol 


7.49 x 10~ 23 


7.49 x 10" 23 


<> m lt m 


2.39 x 10~ 21 


2.40 x 10~ 21 
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Yukawa Coupling Matrix for Point A: 



3.9502 -1.6538 0.5889 
-0.7276 -2.0054 -3.9004 
-1.4061 1.4756 1.5318 



Yukawa Coupling Matrix for Point B: 



3.3847 1.8639 -1.3814 
-1.8107 -0.7219 -0.9499 
-2.5435 -1.0497 -3.3588 



Yl 



N 



-3.8918 -3.9447 -3.8380 

-2.6439 2.5796 3.9962 

-0.9223 -1.3577 0.6417 

2.4435 -1.8006 -1.9575 

0.4198 -3.1594 3.5905 

-0.2505 1.3172 2.1521 



(20) 



(21) 



D 





FIG. 7: The distribution of electron Yukawa couplings (Ye) which give a 'good fit' to the fermion mass data 
in the Dirac case. Neutrinos are assumed to follow normal hierarchy in this analysis. The binning is done 
with an interval of 0.2 



C. Bulk Majorana mass term 

Singlet neutrinos typically accommodate Majorana mass terms in addition to the Dirac mass 
terms. These bare mass terms which break lepton number at a very high scale play an essential 
role in the standard four dimensional seesaw mechanism to generate light neutrino masses. The 
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seesaw mechanism with bulk Majorana mass terms has been first considered in [29J. There have 
been other works which have considered brane localised Majorana mass terms |35 | l4] 4l4"Tll48| . Our 
analysis follows the work of [!?§] and extends it by computing the numerical solutions. The part of 
the action which contains the singlet right handed neutrinos is given by 

S N = j d A x J dy^g~(m M NN c + m D NN + 5(y - irR)Y N LHN) (22) 

where N c = C§N T with C5 being the five-dimensional charge conjugation matrix 10 and itlm = CMk, 
with k being the reduced Planck scale 11 . The bulk Dirac mass for the right handed neutrino is 
parametrized as run = cjyk. As before we consider all the mass parameters to be real. The bulk 
singlet fields N have the following KK expansions: 



00 _. 00 
N L (x,y) = J2-^ w e^NP(x)g^\y) ; N R (x,y) = £ -^e^N^^iy), 

(23) 



n=0 v nR * „ tl V-h> 



2 ^ f_ R dy e° (g^g^ + g^g™) = 6^ (24) 



where gi and gn are profiles of the singlet neutrinos in the bulk. They follow the following 
orthonormal conditions 

iR 
l-irR 

Using this, the eigenvalues equations for the gi,R fields become [29] 

{d y + m D )g { £\y) = m n e° g^' \y) - m M g^\y) 
(-dy + m D )gP(y) = m n e a g < l\y) - m M g { [\y) (25) 

where we have assumed the five dimensional wave functions to be real. Unlike the Dirac and higher 



dimensional LHLH term cases, the present system of equations, in Eq.(25) are not consistent with 
a zero mode solution m n = for mn 7^ 0. This is because the zero mode solutions, oc ev c n~ c m 
do not satisfy either Dirichlet or the more general {dy + m^g^y) = boundary condition. Thus in 
the following analysis, we will consider the first KK mode not to be the zero mode but m n = mny 



Furthermore, Eq.(25) does not have simple analytical solutions, though numerical solutions exist. 



We have obtained the numerical solutions of gi,R by solving the second order equations derived 



from Eq.(25). The equation for the Z2 even part takes the form: 



„ m n kRe kR y , ( c N m n e kR yk 2 2 2 / kR \2\ „ 

9l{v) ~ m n e^ - c M k 9dv) ~ { m n e^y-c M k + ^ ~ ~ J R 9l{v) = ° (26) 



10 C5 is taken to be C4. 

11 Majorana mass terms does not have the same interpretation in the bulk as in 4D. 
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FIG. 8: The distribution of neutrino Yukawa couplings (Y^) which give a 'good fit' to the fermion mass 
data in the Dirac case. Neutrinos are assumed to follow normal hierarchy in this analysis. The binning is 
done with an interval of 0.2 

The second order equation for the Z2 odd part is given as 



g'kiv) 



m 



n kRe kRy 



-Mv) 



-c N m n e 



+ c 2 N k 2 - (m n e kR y - c M k) 2 R 2 g R (y) = 0, (27) 



m n e kR y — cuk \m n e kR y — cuk 

where we have used the notation mo = c^k and Mm = cuk introduced earlier. The primes on 
gL{y) and gn{y) indicate derivatives on the profiles. For a given choice of cn and cm one would 
expect to numerically find solutions using the above equations for gL,R as long as they satisfy two 
conditions: (i) wi(i) is also fixed such that the boundary conditions are satisfied consistently (ii) 
There are no singularities in coefficients of the differential equations in the interval [0,7ri?]. This 
second condition requires that for unique solutions, only those values of cm and m n are allowed for 
which m n e a — tum is non zero. Note that this condition is always true when cm is negative. For 
positive cm, the allowed region is shown in Fig. M, where all the shaded region has m n e u — tum 
non zero. As can be seen from the figure, as cm increases, the KK mass scale also increases. In 



Fig. 10 we show solutions to Eq.(26) for a fixed value of cat = 0.58. cm is varied from 0.55 to 



1. From the figure, it is clear as the profile becomes oscillatory as cm becomes greater than cat. 
In fact the solutions are sinusoidal for ca/=1 and cat=0. We now address the question of fitting 
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FIG. 9: Region of cy[-m\ parameter space, for positive cm for which the coefficients of the differential 



equation in Eq.(26l are analytic in the interval [0, TtR] 



the lepton masses and mixing. The charged lepton mass matrix has the same form as in earlier 
sections 

Y E 



(0,0) v ~ v 



, r,YE = j^ N (c L )N (c E ) e 



(l-c L -c E )kRir 



(28) 



KK 



Choosing g^' to be the Z2 even profile for the right handed neutrino, the Dirac mass matrix takes 
the form 



m 



(0,1) _ Yn_ 
D ~ Rtt 



N (c L )e^- c ^gi\nR) 



(29) 



where g9 (y) is the solution to Eq.(26). The singlet Majorana mass matrix is determined in the 
flavor space by the choice of cat and cm for each of the generations. For simplicity, for the present 
analysis, we take all of them equal cat, = cat and cm^ = cm for all the three generations 12 . With 
this, singlet neutrino mass matrix becomes proportional to the unit matrix Mr = 1 m^y The 
light neutrino mass matrix now takes the see-saw form given by 



1 tm\T 



m (o,o) =m (o,i)^ m (o,i) +Q 



( ( (o,fc) 



V 



rn 



(k) 



(30) 



where higher order corrections are from higher KK states. To fit the neutrino masses and mixing 
angles we neglect higher order corrections as before. Defining = 2/cY/v, we have 



m£ '°> = Y'^-^^gL^R^M^Y^-^^gL^R) 



(31) 



This can be achieved by imposing an 0(3) symmetry on the N fields. 
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FIG. 10: The Figure shows the form of the profile for solution to Eq.(26 1 for a fixed bulk dirac mass of 0.58 



for the right handed neutrinos. We see that the profile becomes oscillatory as cm becomes greater than cj\r. 



In Table (VI), we present two sample points one for inverted hierarchy and another for normal 
hierarchy, which fit the neutrino masses and mixing angles as well as charged lepton masses with the 
accuracy we have specified in section |nj Both these examples 13 have cm < cjv- The corresponding 
Yukawa coupling matrices are presented in Eqs. ( 32|33 ). 
Yukawa parameters for inverted hierarchy 



2.73 1.81 .108 
-0.83 0.975 .328 
0.327 -0.679 .182 



3.44 -0.41 .87 
0.62 1.583 0.332 
2.74 0.55 2.33 



(32) 



These solutions require that the profiles of the N fields have very small values on the UV brane. 
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TABLE VI: Sample points with corresponding fits of observables for Normal and Inverted Hierarchy schemes 
in Bulk Majorana case with 0(1) Yukawas. The masses are in GeV 



Parameter 


Normal 


Inverted 


M kk 


161.4 


161.4 




0.55 


0.55 




3 x 1(T 13 


1.2 x 10~ 12 




0.58 


0.59 




0.56 


0.57 


cl 3 


0.55 


0.55 


ce ± 


0.735 


0.735 


C E 2 


0.5755 


0.575 


CE 3 


0.501 


0.501 


c m 


0.58 


0.58 


m e 


5.09 x 10~ 4 


5.08 x 10~ 4 


"V 


0.1055 


0.1055 


m T 


1.77 


1.774 


#12 


0.58 


0.58 


#23 


0.80 


0.8 


#13 


0.13 


0.13 


a^L 


7.8 x 10~ 23 


7.8 x 10~ 23 




2.4 x 10~ 21 


2.4 x 10~ 21 



Yukawa parameters for normal hierarchy 



2.56 1.69 1.26 
-0.795 0.927 3.89 
0.414 -0.859 2.86 



Y' E 



2.825 -0.41 .87 
0.62 1.2008 0.332 
2.74 0.55 2.31 



(33) 



D. Brane localized Majorana mass term 



Following our discussion with a bulk Majorana mass term, there could be special cases where 
the Majorana mass term could be localized on either boundary. In this case the bulk profiles 
for the right handed singlets iVj remain unchanged. The eigenvalue equations are same as in Eq.©. 
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1. UV localized mass term 



The case with UV localized Majorana mass term was studied in 
case is given as 



. The action in this 



S N = j d A x J dy s[=g (5(y)NN c + m D NN + 5(y - irR)Y N LHN) 



(34) 



where we have expressed ium = S(y). Substituting the KK expansions from Eq.(17), the effective 
4-D neutrino mass matrix, in the basis \ T = { v ^l\ j takes the form 



-X T M NX 



where M 



(0,0) . 



M N 



is defined in Eq.(19). Let f^(0) denote the value of the profile of the first KK mode 






M°' 0) 







a4°' 0) 











M Maj 


M Maj 


M K k 








M KK 






(35) 



of N at the UV brane i.e, y=0 and /jv(0), defined in Eq.Q, is the zero mode profile of N evaluated 

,(o,i) _ 



at y=0. The individual elements of Eq.(35) are then defined as: Ai)?' 1 ' = y=/jv(7T.R)Y^; 
M%o% = tjs&iO) ;M%% = &f$\0)fir(py, M™% = (0)^(0) and M KK is the KK 

mass of first KK mode of N. The small neutrino masses can be fit by choosing cn ~ 0.32 for which 



M 



Maj 
,(0,0) 



10 GeV. The charged leptons are fit by choosing cl,e > 0.5. This scenario along with 



flavour implications has been extensively dealt in [35J. 



2. Pure Majorana Case 

An interesting sub case of the Bulk Majorana term would be the situation where ran = c^k = 0. 
As we have seen from the discussion in the previous section, in such a case, the profile equations 
become oscillatory. The eigenvalue equations now take the form: 

d v9£\y) = m n e a g < £\y) -rriMg^iy) 
-d y gP{y) = m n e CJ g { jl\y)-m M g ( £\y) (36) 

Contrary to the Dirac+ Majorana case of the previous section, the above set of equations allow 
solutions for zero modes, tjiq = 0. The solutions are given as 

9l(v) = Ncosi^Y mMy ) 

9n(y) = A r sin(^^ m M y), (37) 
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where N is the normalization factor given by N = VTrRke~°- 5a( - nR \ These solutions are consistent 
with the boundary conditions. The neutrino mass matrix has a specific structure in this case, 
as there are contributions from the first KK mode, which might be important. In the basis, 
X T = {u^\N^\N^} the mass matrix takes the form 

rrin' ^ m^' 1 ^ ^ 



C m = -^ x T Mx ; M 



mg' 0) 
\m D m (1) J 



(38) 



From the above, we see that at the zeroth level, light neutrino and singlet neutrinos form a pseudo- 
Dirac structure, leading to maximal mixing between these two states. For the three flavor states, 
we would have three light states which are sterile. We have not pursued the phenomenology of this 
model further. 

III. LEPTON FLAVOR VIOLATION 

We now study lepton flavor violating constraints on the three neutrino mass models considered 
in the present work. Lepton flavor violation within the RS framework has been studied in detail in 
|30| . The localization of the fermions in the bulk at different places leads to non-zero flavour mixing 
between the zero mode SM fermions and higher KK states, which contribute to flavor violating 
processes both at the tree and the loop level. The tree level flavor violating decay modes of the 
form li — > Ijlkh are due to non-universal overlap of the zero mode fermions with the Z-boson KK 
modes. At the 1-loop level, penguin graphs contribute to rare decays like lj — )■ li + 7. The SM 
states mix with their heavier KK states on the IR brane, and thus may give rise to significant 
contributions to dipole processes in particular. The present LFV limits are very strong and are 



listed in Table VII 



In this section we calculate the Branching fractions for the leptonic FCNC. The effective 4-D 
lagrangian describing I — > I' process is given by [30] 

+ a 5 (ei ? 7 /i /x R ) {eLl^h) + a 6 (e L ^fi L )(e R 'y^e R )] + h.c. (39) 
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TABLE VII: Present Experimental Bounds on LFV Processes 



Process 


Experiment 


Present upper 


bound 


BR(/j e 7 ) 


MEG EMEU] 


2.4 x 10- 


12 


BR(/i -> eee) 


meg gsisn] 


1.0 x 10" 


12 


L.K(fi -> em li) 


CTl\TT^UTTA/r TT IK i 1 

bllNDKUM-li plj 


0.1 X ID 


13 


BR(r 


BABAR/Belle [5J] 


4.4 x 10" 


-8 


BR(r e 7 ) 


BABAR/Belle [55] 


3.3 x 10" 


-8 


BR(t — > /Lt /x /x) 


BABAR/Belle [5J] 


2.0 x 10" 


-8 


BR(t -^eee) 


BABAR/Belle [5J] 


2.6 x 10" 


-8 



A. Tree level decays 

The breaking of the electroweak symmetry at the IR brane mixes the zero mode gauge boson 
with the higher modes. To parametrize this mixing, let (Z(°), ZW) and (Z /(0) denote the 

gauge boson states before and after diagonalisation of the gauge boson mass matrix respectively. 
Assuming only one KK mode for simplicity, they are related as [30] 

(40) 

where M z (t) is the mass of first KK excitation of the Z boson. Owing to its flat profile the Z^ 1 
couples universally to all three generations. However, the coupling of Z^ l \ whose profile is peaked 
near the IR brane, is generation dependent. This coupling depends on the localization of the 
fermions along the extra-dimension thus giving rise to non- universality. Let rj T = {eu^M-, t m} 
be vector of fermions in the mass basis. Let c£' be a 3 x 3 matrix which denotes the coupling of 
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SM fermions in the mass basis to Z 1 



(i) 



It is given as 

h 



L,R = 9L,R r)L,R-Dl R . 



I„ o 



■Dl,r-Vl,r i' 



(i) 



(41) 



J T 

where <7l,.r is the SM coupling, D^r are 3x3 unitary matrices for rotating the zero mode (SM) 
fermions from the flavour basis to the mass basis. I is the overlap of the profiles of two zero mode 
fermions and first KK gauge boson. It is given by 

He) = £ R dye^(tf\y,c)) 2 ^(y) (42) 

£^'(y) denotes the profile of the first KK gauge boson. It is plotted as a function of a generic 
bulk mass parameter c in Fig.[l3]. As we can see from this figure, the overlap function 1(c) 
becomes universal for c > 0.5 and for c < —15. The off diagonal elements of a$ represent the 
flavour violating couplings. The contribution to li — > Ijhh from direct exchange is suppressed 
compared to that of Z^°\ The contributions to the coefficients a£ 6 in Eq.(39) due to the flavour 
violating coupling of Z^ as well as direct Z^> exchange are given as 



2 /// ,.;v2/,/ 1 - lj)^ , 



(43) 



'z(i) 



-2g L (V2kRn - Ij) 
-2g L (V2kRn - Ij) 



vrir 



,(!)« 



M 2 L 



vrir 



M 2 R 



-2g R (V2kRir - Ij)-^-a 



M 2 " L 



The Branching fractions for the tree level decays are given as |3U| 
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(44) 



Similarly, the relevant quantities for jj, — > e conversion in Ti are given as: 
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(45) 
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FIG. 11: Tree level contribution to /i — > eee due to exchange of Z'^ 1 '. The effective contribution is 
proportional to this graph. 

where p e ~ E e ~ m^. Gf is the Fermi constant, a is the electromagnetic coupling. The most 
stringent constraint for fi — e conversion comes from Titanium (Ti^)- Atomic constants are defined 
as: Qn = v u (2Z + N) + v d (2N + Z) with N being the neutron number, Z e ff = 17.61, form factor 
F p = 0.55, T capt = 2.6 x 10 6 s' 1 for Titanium [S3]. 



B. Dipole Transition lj — » ^7 



The dominant graph is due to scalar exchange in the loop. One of them is due to Higgs exchange 



as shown in Fig. 12 . The amplitude for this process is given as 



AI 



(2vr) 4 tK * J L E p' 2 -M 2 p 2 -M 2 E p 2 -Ml E 



k 2 — m 2 H 11 



(46) 

where p = p — k, p' = p' — k and q = p — p'. A l L R = F % L r Dl^r and M n denotes the mass of the n th 
mode KK fermion. Fl e is a function of bulk masses which are taken to be diagonal in the flavour 
space. It is given as 



Fl.r 



fc Ll ,c Bl {nR) 
/, 








cl 9 ,ce. 



(irR) 






fc L ,,c E AnR) 



(47) 



The amplitude for Eq.(46) can be rewritten as 



M (j -> j 7 ) = {eDlF L Y^vY^F R D R )ijJ(p,if,q) 



(48) 



The expression J(p,p' , q) is the momentum integral in Eq.(46 ). It is log divergent owing to a double 



independent sum over two KK modes. We regularise it using a cutoff of A ~ 4vrM^ } ~ 15 TeV. 
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FIG. 12: Higgs mediated j — ¥ i"f. The dot represents the mass insertion. Flavour indices have been 
suppressed in the internal charged KK lines. (L,R) represents the KK modes corresponding to the left and 
right chiral zero modes. 

The other dominant contribution is due to Fig.[l9] is discussed in Appendix[B]. The Branching 
fraction for the dipole decays lj — > is given as 

1 2"7T 2 

BRili "> Un) = jGrm^^ 1 + A *> (49) 



where the coefficient due to Figs. 12|19 is given as 



a l = 2 ^j4^ FiMt+y ^ t)Wii (5o) 

and Ar = A\. The other dipole contributions are discussed in Appendix |b]. We now proceed to 
discuss the LFV rates for the mass models discussed in Section [nj. The quantities, like the KK 
masses of fermions, the rotation matrices Dl^r etc. which determine the LFV rates are functions 
of the bulk mass parameters. We compute these quantities for each point of the best fit parameter 
space obtained earlier for the LHLH and the Dirac case and use it to constrain the parameter space 
from flavour violation. 



C. LHLH Case 

The contributions to trilepton decays from graphs like Fig. [IT] are highly suppressed in the 
parameter space of interest. This is because the couplings of the zero mode fermions to the KK 
gauge boson become universal for the fermions sufficiently localized towards IR and UV branes, 
as can be seen in Fig.[l3]. However, there could be other potentially large contributions. This 
comes from the large mixing between zero mode charged singlet states and the first KK modes of 
the lepton doublets; the corresponding Yukawa coupling is very large due to the large negative eg 



values. Example of such a graph is shown in Fig. 14 . Exact value of the contribution, of course 
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FIG. 13: Coupling of two zero mode fermions to Z\ as a function of bulk mass parameter [? 



depends on the values of Dl,r and other parameters. We have not considered these graphs in the 
present work. We note that for a fairly degenerate bulk doublet masses, (cxj, the combination of 
the matrices which enter in these graphs are aligned with the zero mode mass matrix for charged 
leptons. The best parameter space does contain such regions where all the are degenerate. 
We found several examples of that kind. Another potential problem with the highly localized 
IR charged singlets, is the shift in the universal coupling constant g^. This could effect Z —> 11 
branching fractions. Models with custodial symmetries or very heavy KK gauge bosons could avoid 
this problem. We have not addressed this issue here. 
Finally, contribution to U 
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are heavily 



Zj7 due to loop diagrams of the form in Fig. 
suppressed owing to the heavy KK mass scales corresponding to the charged singlets. The corre- 
sponding masses are in shown in Table |TT]. Additionally, the large effective 4-D Yukawa couplings of 
the charged singlets to the KK modes make it difficult to apply techniques of perturbation theory 
to calculate graphs like those in Fig. 12|19 . 



D. Constraints on Dirac Neutrinos 



The Dirac case gives a good fit to the leptonic data for a reasonable choice of 0(1) parameters. 
However, the parameter space is strongly constrained from flavour considerations. In the parameter 



space of interest the dominant contribution to tree-level decays comes from Fig. 11 . The parameter 



space of the bulk doublets and charged singlets consistent with tree level contribution is shown 



in Fig. 15 . The lightest M Z {X) mass required to satisfy all constraints from tree-level processes 
~ 1.9 TeV. Fig.|15| shows the points within the best fit parameter space consistent with all 
constraints from tree-level processes. As can be seen from the figure, very few points pass the 



31 




FIG. 14: Additional tree level contribution to /.i — > eee. For a fairly degenerate bulk doublet mass in LHLH 
case this contribution is negligible. For the Dirac case this graph receives wave function suppression in 
addition to the KK scale suppression. 

constraints. The black point is allowed for a KK gauge boson scale of 1.9 TeV, where as the 
green points are for mass of 3 TeV. The constraints from dipole processes are far more severe. 
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FIG. 15: The black dot and the green region represent the parameter space permitted by tree-level constraints 
for a KK gauge boson scale of 1920 and 3000 GeV respectively 



Corresponding to the ex e values in the best fit parameter space, the mass of the first KK excitation 
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of the leptons varied from approximately 850 GeV to 1400 GeV as presented in Table (IV). We 
found no points which satisfied the constraints from /i — > ej, r — > ej and r — > simultaneously. 
The constraint from /i — > ej was most severe and required a KK fermion mass scale O(10) TeV to 
suppress it below the experimental limit given in Table (VTl] . 



E. Constraints on scenarios with bulk Majorana mass 

The tree-level decays only constrain the parameter space of the bulk doublets and charged 



singlets as seen in Fig. 15 . Since, the charged lepton mass fitting is independent of any right 
handed neutrino parameter, the constraints coming from tree-level decays in the Dirac case are 
applicable in this case as well. 

The contribution to dipole decays of the form lj — > lij due to charged Higgs shown in Fig. [19] 

g^^R) is required to be small to fit neutrino 



is small. This is because, as shown in Table 



VI 



masses. Thus, the dominant contribution to dipole decays in this case is due to Higgs exchange 
diagram shown in Fig. fl2| . They are calculated for the both the normal and inverted hierarchy 



cases presented earlier and are given in Table VIII . The branching fractions are evaluated for 



Mkk ~ 1250 GeV which is the first KK scale of the doublet. 

TABLE VIII: BR for dipole decays for the case with bulk Majorana mass 



Hierarchy 


BR(^ ->■ ej) 


BR(r /i7) 


BR(r -> e 7 ) 


Inverted 
Normal 


2.4 x 1CT 5 
1.4 x 1(T 5 


1.9 x 1CT 5 
3.4 x 1CT 5 


7.6 x 10" 6 
1.3 x 1CT 5 



IV. MINIMAL FLAVOR VIOLATION(MFV) 

From the discussion above it is clear that lepton flavor violating constraints are strong on RS 
models with fermions localized in bulk and Higgs localized on the IR brane. In the Dirac and the 
Bulk Seesaw case flavor violation rules out most of the 'best fit' parameter space. One option to 
evade these bounds would be to increase the scale of KK masses. As we have seen in the LHLH 
case, the fits indicate to the highly hierarchal spectrum with KK masses of the O(10 2 ) TeV for the 
singlet charged leptons, the flavor violating amplitudes are highly suppressed and thus do not put 
severe constraint on the model. However, the Dirac and the Majorana cases whose best fit regions 
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have lighter KK spectrum would essentially be ruled out. The misalignment between the Yukawa 
coupling matrix and bulk mass terms which determine the profile is the cause of the large flavor 
violating transitions leading to strong restrictions on these models. In |54j the authors imposed 
discrete symmetries to constrain Flavour Changing Neutral Currents (FCNC). In this work we 
adopt the Minimal Flavour violation ansatz which reduces the misalignment by demanding an 
alignment between the Yukawa matrices and the bulk parameters. 

The ansatz of Minimal Flavour violation was first proposed for the hadronic sector |55| . It 
proposes that new physics adds no new flavor structures and thus entire flavor structure in Nature 
is determined by the Standard Model Yukawa couplings. In the leptonic sector, MFV in not 
uniquely defined due to the possibility of the seesaw mechanism. Several schemes of leptonic 
minimal flavor violation are possible [56J. 

The proposal to use the MFV hypothesis in RS was first introduced in |31] in the quark sector. 
There were subsequent extensions in the leptonic sector by |32l [35] . The MFV ansatz assumes that 
the Yukawa couplings are the only sources of flavor violation. In the RS setting this would require 
that the bulk mass terms should now be expressed in terms of the Yukawa couplings [31] . The 
exact expression would depend on the particle content and the flavor symmetry assumed. 

A. Dirac Neutrino Case 

In the presence of right handed neutrinos the flavour group is SU(3)l x SU(3)e x SU(3)n', 
the lepton number is conserved. The Ye transforms as Ye —> (3, 3, 1) and Y/v transforms as 
Yat —> (3,1,3). The Yukawa couplings are aligned with the five dimensional bulk mass matrices. 
The bulk masses can be expressed in terms of the Yukawas as 

c L = atl + a 2 Y' ' e y£ + a 3 Y^Y$ c E = bY^Y E c N = cY$Y' N (51) 

where a,b,c G 5i and Y' E N are as defined earlier as Y' E N = 2kY E ^. Owing to the flavor symmetry 
we work in a basis in which Y' E is diagonal. We then rotate Y^ by the PMNS matrix i.e, writing 
Yn -> Vp MN sBi&g(Y^) where the Diag(Y^) = Diag(0. 709, 0.709, 0.75). The c L value chosen 
is 0.5802 for all three generations. The cjv values chosen are respectively 1.17241, 1.172, 1.311 
respectively. The bulk singlet mass parameters are eg = (0.7477,0.58059,0.401) 

The simplest Yukawa combination transforming as (8,1,1) under the flavour group is given as 

A = Y^Y* (52) 
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Thus the BR for fj, — )■ cy, which is the most constrained is given as |35j 

BRfa 4e 7 )=4x 10- 8 {Y> N Y$)\ 2 (53) 



Y'n 



0.586033 0.383951 0.115044 
-0.335962 0.370429 0.53165 
0.215349 -0.466953 0.516346 



(54) 



The (1,2) element of A which is responsible for — > e^f is 0.006 which gives a contribution of 
1.44 x 10 -12 , for a fermion KK mass of around 3 TeV. 



B. Bulk Majorana mass term 



Owing to the presence of a bulk Majorana mass term, we choose the flavour group for the 
lagrangian in Eq.Q is SU(3) L x SU(3) E * 0(3) N . Y E transforms as Y E -»• (3,3,1) and Y N 
transforms as Y/v ~~ > (3,1,3). The bulk Majorana term N C N transforms as (1,1,6) under this 
flavour group. In terms of the dimensionless Yukawa couplings, Y' E N the bulk mass parameters 
can be expressed as 

c L = a 1 I + a 2 Y , E Y* + a 3 Y^Y{f c E = 1 + bY$Y E c n = 1 + cY^Y}, c M = dl 3x3 (55) 

where a,b,c,d G K. cm = 0.55 and cn = 0.58 are chosen for the right handed neutrino bulk 
mass parameters. The value of profiles for the singlets are chosen appropriately at the boundary so 
as to fit the neutrino data using the 0(1) Yukawa couplings. As before we work in a basis in which 
Y' E is diagonal. In this basis Y' N = VpMATS'Diag(Y^). This removes the dominant contribution to 



dipole decays due to the Higgs exchange in Fig. [12] . The contribution due to Fig. 19 is very small 
owing to wavefunction suppression of the singlet neutrinos. Thus, we see that the MFV ansatz is 
successful in suppressing FCNC's for both the Dirac and the bulk Majorana case. 



V. SUMMARY AND OUTLOOK 



Understanding neutrino masses and mixing is an important aspect of most physics beyond the 
Standard Model frameworks. The Randall-Sundrum setup while solving the hierarchy problem 
could also form a natural setting to explain flavour structure of the Standard Model Yukawa 
couplings. The quark sector has already been explored in this context in detail. While there have 
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been several analysis in the leptonic sector, in the present work we have tried to explore the same 
in a comprehensive manner, filling the gaps wherever we found it necessary. Our aim had been to 
determine quantitavely the parameter space of both the 0(1) (dimensionless) Yukawa couplings as 
well as the bulk mass parameters which can give good fits to the leptonic data. 

We have concentrated on the RS setup with the Higgs field localized on the IR boundary. We 
have considered three cases of neutrino mass models (a) The LH LH higher dimensional operator 
(b) The Dirac case and the (c) Majorana case. The LHLH fits require large negative c-parameters 
which reflect the composite nature of the charged singlets. There is some parameter space in this 
case where the flavor constraints are weak. However, the model has very large effective 4-D Yukawa 
couplings between the zero mode SM fermions and the KK fermions, which makes it unattractive 
from perturbation theory point of view. We have also presented the distributions of the Yukawa 
couplings in the best fit region. Most of the individual Yukawa couplings are concentrated on 
the higher side of the 0(1) range we have chosen. The Dirac and Majorana cases offer large 
parameter space without the need of large hierarchies in the c parameters. We have also presented 
the distribution of the Yukawa couplings in the Dirac case. We could not find strong correlations 
between the Yukawa couplings and the c-parameters. There are strong constraints from the lepton 
flavor violating rare processes. These can be circumvented by a suitable choice of Yukawa couplings 
and c-parameters guided by the MFV ansatz. The Majorana case, in particular allows for several 
classes of MFV schemes, which will be explored in an upcoming publication |57j . 

While we restricted ourselves to the Higgs located on the IR brane, it can also be allowed to 
propagate in the bulk. Lepton flavor violating amplitudes however are now cut-off independent, 
which makes the computations more predictive. But with the Higgs boson in the bulk one has to 
invoke other scenarios like supersymmetry to solve the hierarchy problem. 
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APPENDIX A: INVERTED MASS FITS 



We present the results of the scan performed for inverted hierarchy for both the LLHH and 
the Dirac case. In the case for the normal hierarchy it was easier to find c values and order 
one Yukawa entries which satisfied all constraints. However, the choice of these parameters 
which fits the data in the inverted case is very subtle. This is because one requires two large 
mass eigenvalues in the inverted case which must satisfy the Am^ constraint. This requires 
a very careful choice of order one Yukawa parameters. The parameter space for c values does 
not differ much between the normal and the inverted case. For the case of inverted hierarchy, 
we choose points which satisfy < % 2 < 10. For the Dirac case we performed a scan only for c > 0.5. 

(A) LHLH case 

TABLE IX: Sample points for Inverted Hierarchy in LHLH case with 0(1) Yukawas. The masses are in GeV 



Point 


A 


B 


x 2 


7.48 


6.61 


CLi 


0.8967 


0.9162 


CL 2 


0.8983 


0.8920 


CL 3 


0.8913 


0.8945 


CEi 


-3758.1502 


-2099.8993 


ce 2 


-6005847.4955 


-552577.8188 


ce 3 


-32730342.0982 


-23953472.2265 


m e 


5.11 x 10~ 4 


5.09 x 10~ 4 




0.1056 


0.1056 


m T 


1.775 


1.755 


$12 


0.584 


0.55 


$23 


0.829 


0.875 


$13 


0.148 


0.160 




7.49 x 10~ 23 


7.46 x 10~ 23 


Sm 2 atm 


1.90 x 10~ 21 


2.7 x 10~ 21 
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Yukawa for Point A 



0.8249 0.8516 1.1111 
1.3600 1.5956 1.8402 
3.5831 3.5664 2.9092 



Yukawa for Point B 



Y' E 



2.5874 0.5123 3.6064 
3.9696 2.4876 1.9903 
3.8604 1.1438 3.9712 



(B) Dirac Case 



-3.5528 2.6612 1.4503 
2.6612 3.8149 1.2903 
1.4503 1.2903 -0.6682 



(Al) 



-3.6860 -3.6778 3.9987 
-3.6778 2.1362 3.3252 
3.9987 3.3252 -0.8497 



(A2) 



TABLE X: Sample points for Inverted Hierarchy in Dirac case with 0(1) Yukawas. The masses are in GeV 



Parameter 


Point A 


Point B 


x 2 


0.30 


8.04 




0.5565 


0.51 




0.5556 


0.5316 


c La 


0.5433 


0.5012 


CEi 


0.7681 


0.8092 


ce 2 


0.6186 


0.6498 


ce 3 


0.5044 


0.5674 




1.2450 


1.2765 


c N2 


1.2421 


1.2755 


c N3 


1.2546 


1.2941 


m e 


5.1 x 10~ 4 


5.08 x 10~ 4 


«v 


0.1055 


0.1055 


m T 


1.769 


1.81 


#12 


0.59 


0.59 


#23 


0.80 


0.72 


#13 


0.155 


0.152 




7.49 x 10~ 23 


7.48 x 10~ 23 


Sm 2 atm 


2.40 x 10~ 21 


2.16 x 10~ 21 
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Yukawa for Point A 



Y' 



2.2645 2.7691 0.4272 
1.0499 -3.6695 -1.0818 
-2.2402 -0.5400 -1.9176 



Yukawa for Point B 



-3.7916 -0.3960 -2.5573 
1.2699 -2.3757 3.2167 
-3.5010 3.4430 2.8224 



Y' 



N 



-0.2202 -2.3054 1.5602 
3.4794 -2.2140 0.2302 
-2.0676 -1.7529 0.7888 

-3.9443 -0.9714 0.1848 
-2.5788 0.2609 3.3684 
0.5020 -3.0268 -3.1765 



(A3) 



(A4) 
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FIG. 16: The plot represents the parameter space for the bulk masses of charged doublets for inverted 
hierarchy 
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FIG. 17: The plot represents the parameter space for the bulk masses of charged singlets for inverted 
hierarchy 
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FIG. 18: The plot represents the parameter space for the bulk masses of neutrino singlets for inverted 
hierarchy 

APPENDIX B: AMPLITUDES FOR DIPOLE TRANSITIONS 

In this section we review the other potential contributions to the dipole processes % — > jj 
a) Internal flip in neutrino KK line in Dirac Case 



This contribution is absent for the LHLH case as it involves neutral internal KK lines cor- 
responding to the right handed neutrino. In the unitary gauge the charged Higgs is nothing but 



the longitudnal component of the W boson. This displays a similar divergence to Fig. 12 owing 




TO 



FIG. 19: "Charged" Higgs mediated j — > ij. The dot represents the mass insertion. Flavour indices have 
been suppressed in the internal neutral KK lines. (L,R) represents the KK modes corresponding to the left 
and right chiral zero modes respectively. 
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to the presence of double KK sum. 



(wv , v ,i v , r d 4 k _ , ^ {P' + M n ) p + M n 1 1 

= (F L Y N Y N e T Y E F E ) l3 J £ —^(p )(2F-<f ) ^ _ ^ ^ _ ^ _ ^ (fc _ g)2 _ m y 

(Bl) 



b) Gauge contribution 




FIG. 20: Contribution to the dipole graph due exchange of KK gauge bosons and charged KK fermion lines. 



The amplitude for Fig. 20 



is given as 



M . . - ( A Q,n,l v Y > A Q, m ,l, f dk I + M n p + M n p' + M m 



ii.ru 



(B2) 

where A '™' 1 represents the coupling of zero mode fermion to n th mode fermion and I th mode 
gauge boson. The contribution from this sector is suppressed in both the Dirac and LHLH case in 
the parameter space under consideration. 
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